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The presence of sulfur atoms on the Pd(100) surface is 
known to hinder the dissociative adsorption of hydrogen. Us- 
ing density-functional theory and the full-potential linear aug- 
mented plane-wave method, we investigate the potential en- 
ergy surface (PES) of the dissociative adsorption of H2 on 
the sulfur covered Pd(100) surface. The PES is changed sig- 
nificantly compared to the dissociation on the clean Pd(100) 
surface, in particular for hydrogen close to the S atoms. While 
the hydrogen dissociation at the clean Pd(100) surface is non- 
activated, for the (2x2) sulfur adlayer (coverage &s= 0.25) 
the dissociation of H2 is inhibited by energy barriers. Their 
heights strongly depend on the distance between the hydro- 
gen and sulfur atoms leading to a highly corrugated PES. The 
largest barriers are in the vicinity of the sulfur atoms due to 
the strong repulsion between sulfur and hydrogen. Still the 
hydrogen dissociation on the (2x2) sulfur covered Pd(100) 
surface is exothermic. Thus the poisoning effect of sulfur 
adatoms for H2 dissociation at low sulfur coverage (Os < 
0.25) is mainly governed by the formation of energy barri- 
ers, not by blocking of the adsorption sites. For the c(2x2) 
sulfur adlayer (Gs= 0.5), the PES for hydrogen dissociation 
is purely repulsive. This is due to the fact that for all dif- 
ferent possible adsorption geometries the hydrogen molecules 
come too close to the sulfur adatoms before the dissociation 
is completed. 

PACS numbers: 68.45. Da, 73.20.At, 82.65.Jv 



I. INTRODUCTION 

The modification of the chemical reactivity and selec- 
tivity of surfaces by adsorbates is important for, e.g., 
building better catalysts. Therefore the investigation of 
the microscopic mechanisms of how adatoms poison or 
promote certain reactions is - besides its fundamental 
interest - of great technological relevance. Hydrogen dis- 
sociation has become the benchmark system for theoreti- 
cal and experimental studies of simple chemical reactions 
on surfaces. On clean metal surfaces, a detailed picture 
of how H-H bonds are broken and how new bonds be- 
tween the hydrogen atoms and the surface-.are formed 
has been developed by theoretical studies.Erla The modi- 
fication of the potential energy surface (PES) of hydrogen 
dissociation on an adlayer-covered surface has not been 
addressed in similiar detail by theory so far. 

Experimentally it is well established that the presence 
of sulfur on metal surfaces causes_a_drastic reduction 
in hydrogen sticking probabilities .lijliia Earlier theoret- 



ical studies explaining the poisoning mechanis m . fo cused 
on general concepts. Feibelman and Hamannotj sug- 
gested that the poisoning effect of sulfur is related to the 
sulfur induced change of the density of states (DOS) at 
the Fermi level. This explanation was also employed by 
MacLaren and co-workers£j Current studies, however, 
have emphasized that the reactivity of surfaces capapt 
be solely understood by the DOS at the Fermi level£Mi!j 
A different mpdel_for the poisoning was proposed by 
N0rskov et al.ilB~c3 These authors explain the modifica- 
tion of the reactivity by adlayers by the interaction of the 
H 2 molecule with the adlayer induced electrostatic field. 
These models, however, do not provide a detailed micro- 
scopic picture of how the hydrogen dissociation process is 
actually modified by the adlayer, i.e., either by blocking 
adsorption sites for atomic hydrogen or by building up 
energy barriers along the dissociation pathways of H2 . 

Recently the poisoning mechanism of Pd(100) by a sul- 
fur adlayer has,b«en investigated theoretically by Wilke 
and Scheffler.E3£3 The Pd(100) surface is a well-suited 
substrat for the investigation of the poisoning mecha- 
nism of catalytic reactions, because many experimen- 
tal and theoretical studies exist and form a wealth 
of information for comparison. At the clean surface, 
hydrogen molecules dissociate spontaneously, i.e., non- 
activated dissociation pathways exist with no hampering 
energy barrier BE3ll3'Ot3 When the surface is covered 
with sulfur, .the H2 sticking probability is significantly 
reduced. E3li3 These experiments show that with increas- 
ing sulfur coverages Os the initial sticking coefficient 
of H2 strongly decreases. This is-true in particular for 
molecules with low kinetic energytLHl (< 0.1 eV) that can 
adsorb only if nonactivated dissociation pathways exist 
and are accessible. For Og « 0.25 the initial sticking co- 
efficient of molecules with low energies is approximately 
two orders o f ma gnitude smaller than the one for the 
clean surface Ji3ii3 indicating that for this sulfur coverage 
non-activated dissociation is nearly completely hindered. 
When the sulfur coverage is increased, the initial sticking 
coefficient reduces even further; at Os ~ 0.5 it is about 
three orders of smaller than at the clean surface £j In ad- 
dition to the significant, decrease of the initial sticking co- 
efficient, TPD studiesEJ observed a decrease of the hydro- 
gen saturation coverage with increasing ©3. Burke and 
MadixEJ therefore concluded that sulfur adatoms sub- 
stantially reduce the hydrogen adsorption energy at sites 
in their vicinity, making these positions unstable against 
associative desorption. An earlier permeation study of 
Comsa, David, and SchumacheiO led to the same con- 
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elusion. 

These findings-are at variance with the results of Wilke 
and Scheffler£3i23 They showed by density-functional 
theory calculations that the hydrogen dissociation at the 
(2x2) sulfur covered Pd(100) surface is still exothermic 
and that the blocking of adsorption sites is therefore only 
of minor importance at low coverages. Adsorbed sulfur 
builds up energy barriers in the entrance channel and 
thus hinders the dissociation. Their results gave a clear 
explanation for the poisoning mechanism. We have now 
extended their work. We have investigated in detail how 
the poisoning of the hydrogen dissociation on Pd(100) by 
the presence of sulfur depends on the position and orien- 
tation of the molecule. Furthermore, we have addressed 
the influence of the sulfur coverage on the poisoning. 

In this paper, using density-functional theory and the 
full-potential linear augmented plane wave (FLAPW) 
method, we calculate the PES of H2 dissociative ad- 
sorption over the (2x2) and the c(2x2) sulfur covered 
Pd(100) surface at different adsorption sites and differ- 
ent orientations, and thus provide a complete picture of 
the poisoning mechanism caused by the adsorbed sulfur 
on Pd(100) surface. .-We, fully confirm the conclusions of 
Wilke and SchefflerE3e3: the poisoning effect of sulfur 
adatoms for H2 dissociation over Pd(100) surface at low 
sulfur coverage (0s < 0.25) is governed by the forma- 
tion of energy barriers. The height of the energy barriers 
depends strongly on the distance between H 2 and the S 
atoms. For high sulfur coverage (0s= 0.5), the PES for 
hydrogen dissociation becomes purely repulsive due to 
the fact that for all possible different adsorption geome- 
tries the hydrogen molecules come too close to the sulfur 
adatoms before the dissociation is completed. The results 
are now sufficiently complete for a six-dimensional quan- 
tum dynamical simulationc3. This requires the represen- 
tation of the ab initio PES by an analytical functionEZI in 
order to interpolate the actually calculated points and to 
achieve a closed expression of the PES. This will be also 
presented in this paper. 

The structure of the paper is as follows. In Section 
II we describe the theoretical method and computational 
details. The potential energy surface for the hydrogen 
dissociation on the (2x2) sulfur covered Pd(100) surface 
at different adsorption sites is presented and discussed 
in Section III. Section IV reports the results of the H 2 
interaction with the c(2x2) sulfur covered Pd(100) sur- 
face. The dependence of the PES on the orientation of 
the molecule is analysed in Section V. After a brief de- 
scription of the analytical representation of the PES suit- 
able for a quantum dynamical simulation in Section VI, 
the paper concludes with a summary in Section VII. 

II. COMPUTATIONAL DETAILS 

We have used density-functional theory together with 
the generalized gradient approximation (GGA)l2a for the 



exchange-correlation functional. The full-potential lin- 
ear augmented plane wave (FLAPW) methodErEj is em- 
ployed for solving the non-relativistic Kohn-Sham equa- 
tions. The FLAPW wave functions in the interstitial 
region are represented using a plane wave expansion up 
to -E C ut = 11 Ry For the potential representation plane 
waves up to E cut — 169 Ry are taken into account due to 
a small muffin-tin radius around the H atoms (rj^ T = 0.37 
A). Inside the muffin-tin spheres the wave functions are 
expanded in spherical harmonics with l max = 10, and 
non-spherical components of the density and potential 
are included up to l m ax = 4. For the k-space integration, 
a grid of 5x5 (or 6x6) uniformly spaced points are used 
in the two-dimensional Brillouin zone of the (2x2) [or 
the c(2x2)] surface unit cell, and we find that for a finer 
mesh the adsorption energies change by at most 30 meV 
for the (2x2) hollow site adsorption geometry. 

We used a supercell geometry and modeled the metal 
surface by three-layer slabs which were separated by a 
12 A vacuum region. Hydrogen and sulfur atoms were 
placed at both sides of the slab. For the geometry of the 
sulfur covered Pd(100) surface, the results reported in 
Ref. |2^ were used. Tests were made for a five-layer slab 
at different adsorption geometries, and the average total 
energy difference from a three-layer slab calculation was 
found to be less than 40 meV (see Table Q). rEhis agrees 
with results reported by Wilke and SchefflerEil who also 
studied the influence of the slab thickness on the poten- 
tial energy for slabs with up to seven layers. They found 
energy differences of less than 30 meV for the barrier 
heights and the adsorption energies at the hollow sites of 
the surface. 

The substrate geometry was kept fixed for the H 2 ad- 
sorption pathways studied. This is appropriate and plau- 
sible because the mass mismatch of H with the substrate 
atoms is significant. The energy zero is taken as the en- 
ergy of the geometry where the molecule is sufficiently far- 
away from the surface (Z = 4.03 A) such that there is 
practically no interaction between the molecule and the 
surface. Zero-point corrections are not included in the 
PES. 



III. H 2 AT THE (2x2) SULFUR COVERED 
Pd(lOO) SURFACE 

A. Determination of the potential energy surface 
(PES) 

Neglecting surface relaxation effects, the potential en- 
ergy surface for the dissociative adsorption of a hydro- 
gen molecule over a sulfur covered Pd(100) surface is six- 
dimensional corresponding to the molecular degrees of 
freedom. The coordinates used in this work are defined 
in Fig. 1. The H2 center of mass position is given by 
three Cartesian coordinates (X, Y, Z) for which the origin 
is chosen as the hollow position of the topmost Pd layer. 
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The two rotational degrees of freedom are given by the 
angle 9 of the molecular axis with the Z-axis (cartwheel 
rotation) and the angle with the X-axis (helicopter 
rotation). The distance between the hydrogen atoms is 
denoted by dfj— h- 

To map this high-dimensional energy surface the 
common strategy is to compute two-dimensional cuts 
through the energy surface, so-called elbow-plots, where 
(X, Y, 9, <p) are fixed, and only c?h-h (the bond length 
of the hydrogen molecule) and Z (its height above the 
surface) are varied. Figure 2 shows the surface unit cell 
for the (2x2) sulfur covered Pd(100) surface. Inside the 
unit cell, there are two bridge sites bi and b2, two hollow 
sites hi and h 2 , the top site above the Pd atom tpd, and 
the top site above the S atom tg- In order to obtain a 
comprehensive information about the adsorption behav- 
iorof H2, we have evaluated the elbow-plots at different 
cuts through the six-dimensional configuration space of 
H 2 where the cuts are defined by the site (X, Y) and the 
molecular orientation (9, (f>). 

B. Results 

To distinguish between the different molecular geome- 
tries, we characterize them by the set of fixed coordinates 
(X, Y, 9, (j>). For example, if one defines a as the unit 
length of the (2x2) unit cell, then (0.5a, 0.5a, 0°, 0°) 
refers to a geometry with the hydrogen molecule placed at 
the hi site in an upright position (molecular axis perpen- 
dicular to the surface); and (0.5a, 0.5a, 90°, cf)) refers to 
a geometry with the molecule at the hi site with its axis 
parallel to surface. The geometry in which the molecule 
is parallel to the surface can also be characterized by the 
positions where the hydrogen atoms will be adsorbed, 
and by that over which the center of mass of the molecule 
is situated. For example (0.5a, 0.25a, 90°, 90°) can be 
denoted by hi-bi-h 2 which refers to a geometry with 
the center of mass over the bridge (bi) position and the 
atoms oriented towards the two hollow sites hi and h 2 
in which they will finally be adsorbed. 

1. 2-D cuts through the PES at hollow sites 

For the PES of the hydrogen dissociation over the 
(2x2) S covered Pd(100) surface (9 S = 0.25), two ad- 
sorption geometries bi-hi-bi and-,bi-h 2 -bi have been 
reported by Wilke and Scheffler.Eil Their results show 
that the presence of a the (2x2) S adlayer on Pd(100) 
changes the dissociation process significantly from be- 
ing non-activated at the clean Pd(100) surface to acti- 
vated dissociation. In their calculations, the dissociation 
pathways corresponding to the bi-hi-bi and bi-h 2 -bi 
geometries at the (2x2) sulfur covered Pd(100) surface 
had energy barriers in the entrance channel of 0.1 and 
0.6 eV, respectively. 



Figure 3(a) shows the re-calculated PES for the bi- 
hi-bi adsorption geometry. The dissociation pathway 
has an energy barrier of 0.11 eV in the entrance-chan- 
nel which is consistent with previous calculations E3 This 
configuration corresponds to the minumum barrier dis- 
sociation pathway. We will show below that the energy 
barriers increase as a function of decreasing distance of 
the hydrogen molecule to the sulfur atoms on the surface. 
For that reason it is not surprising that the bi-hi-bi ad- 
sorption pathway exhibits the lowest energy barrier be- 
cause the hi hollow site is the site on the surface furthest 
away from the sulfur atoms. 

Figure 3(b) shows the PES of the tpd-hi-tpd path- 
way. At the clean Pd(100) surface, this_path has a local 
minimum potential energy of -0.23 eV.Eil The calculated 
results show that upon the adsorption of a (2x2) S ad- 
layer an activation energy barrier towards dissociative 
adsorption of 0.13 eV builds up in the entrance channel 
at (Z, c\h-h) = (1-7 A, 0.78 A), and the minimum at (Z, 
) = (0.98 A, 0.94 A) becomes very shallow with an 
energy of 0.01 eV. For this geometry, the two hydrogen 
atoms are about 3.3 A away from the adsorbed sulfur 
atoms, hence one can expect that there is no direct in- 
teraction between H and S atoms. In such a situation the 
energy increase is due to the sulfur-induced modification 
of the local electronic structure at the surface Pd atomcj, 
as will be confirmed below. 

Figure 3(c) shows the PES for the b 2 -h 2 -b 2 pathway. 
In this geometry the dissociation process becomes purely 
repulsive; it differs from the behavior at the bi-h 2 -bi ge- 
ometry where activated dissociation with an energy bar- 
rier of 0.6 eV is possible.CJ This is easy to understand 
because the hydrogen molecule dissociates towards the 
sulfur atoms at the b 2 -h 2 -b 2 geometry, and the direct 
interaction between H and S atoms causes the dissocia- 
tion process to become purely repulsive. 

Figure 3(d) shows the PES for a geometry with the 
hydrogen molecule adsorbed at the hi site in an upright 
position, i.e. with its molecular axis perpendicular to the 
surface. At the clean surface, there exists a shallow local 
minimum of -0.03 eV at (Z, c2 H -h)= (0.60 A, 0.96 A).H 
As is evident from Fig. 3(c), this pathway becomes re- 
pulsive upon the adsorption of the (2x2) S adlayer and 
the energy at (Z, c?h-h)= (0.60A, 0.96 A) has raised to 
0.28 eV from -0.03 eV. There is still a shallow well, how- 
ever. For this geometry, the two hydrogen atoms directly 
above the hi site are at least 4 A away from the adsorbed 
sulfur atoms, so there is no direct interaction between H 
and S atoms. Thus the increase in the potential energy 
has to be due to the sulfur-induced modification of the 
local electronic structure at this site. 



2. 2-D cuts through the PES at bridge sites 

In addition to the dissociation over the hi and h 2 hol- 
low sites, we also considered geometries where the center 
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of mass of molecule is situated at the bridge site and 
the top sites over the Pd or S atoms. In Figures 4(a) and 
4(b), the PES for the hi-bi-h2 and tpd-bi-tpd pathways 
are presented. 

At the clean surface, the dissociation pathway for hi- 
bi-h2 is non-activated. In Fig. 4(a), the PES for hi-bi- 
h2 geometry indicates that the adsorption of H2 becomes 
an activated process with an energy barrier of 0.16 eV in 
the entrance channel (Z=1.8 A and c?h-h= 0.78 A). This 
energy barrier also results-,from the modification of the 
local electronic structure.^ There exists a second barrier 
of 0.13 eV in the exit channel at (Z, ^h-h )=(0.97 A, 
1.05 A). We believe that this barrier might be due to the 
compensation between the direct attractive interaction of 
the Pd atoms with one H atom and the direct repulsive 
interaction of the S atoms with the other H atom, since 
the top of the energy barrier is situated at the value of 
Z close the sulfur adsorption height (1.24 A). 

Figure 4(b) shows the PES for the tpd-bi-tpd geome- 
try. At the clean surface this dissociation pathway has a 
local energy minimum of -0.2 eV before becoming repul- 
sive. Our calculations show that, upon adsorption of the 
S adlayer, the dissociation process becomes purely repul- 
sive and is different from the behavior at the hi-bi-h2 
geometry where dissociative adsorption over an energy 
barrier of 0.16 eV is possible. Again, this can be under- 
stood by the fact that the hydrogen molecule dissociates 
towards the direction pointing to the Pd atoms and close 
to the S atoms. Thus the direct interaction between the 
H and S atoms leads to the disappearance of the shallow 
minimum along this path (E a d= 0.20 eV) at the clean 
Pd(100) surface. 

3. 2-D cuts through the PES at top sites above Pd or S 
atoms 

The PES for dissociative adsorption over the tpd, tg 
sites near the adsorbed sulfur is significantly different 
from the PES at the hi, Ii2, and bi sites that are further 
away from the adsorbed sulfur. Figure 5(a) shows the 
PES for the h 2 -tpd-h 2 geometry. At the clean surface, 
this dissociation pathway has a local minimum of -0.24 
eV in the entrance channel and an energy barrier of 0.15 
eV in the exit channel [see Fig. 2(c) of Ref. 8]. The ad- 
sorbed sulfur atom is about 2 A away from the adsorption 
site tpd and 1.24 A above the topmost Pd layer. One can 
expect that the hydrogen molecule will interact directly 
the S atom before it reaches the topmost Pd layer; this 
therefore raises the energy of the PES. 

Indeed one can see from Fig. 5(a) that, upon the ad- 
sorption of the S adlayer, the local minimum of PES in 
the entrance channel has disappeared and the energy bar- 
rier in the exit channel has raised to 1.28 eV even though 
the dissociation of hydrogen atoms does not point to- 
wards the adsorbed sulfur atoms. 

Figure 5(b) shows the PES at the ts site within b2-tg- 
b2 geometry. In this geometry, the hydrogen molecule 



directly approaches the S atom on the surface which re- 
sults in a strong repulsion. Seen from Fig. 5(b), we 
found that the PES has increased to more than 1.5 eV 
even when the molecule is still 3 A above the topmost Pd 
layer. However, we found an energy minimum of 1.58 eV 
at (Z, c2h-h)= (2.20 A, 1.90 A). Judging from the fact 
that the hydrogen atoms are 2.20 A above the topmost 
Pd layer, we believe that this energy minimum is com- 
pletely caused by the local bonding properties between 
the H and S atoms. 

From the facts found in Figure 5, one can expect 
that if the projected distance d\\ (parallel to surface) be- 
tween the adsorption site and the adsorbed sulfur atom 
is smaller than 2 A, the dissociation behavior of the hy- 
drogen molecule is dominated by the strong repulsive in- 
teraction between H and S atoms. This can be further 
proved by studying the PES of a frozen H2 molecule (i.e., 
e?H-H=const) above the sulfur covered surface. Figure 6 
shows the PES of a hydrogen molecule with dn-H= 0.76 
A with its center of mass moving inside the (010) plane 
crossing the hi, tpd, ts adsorption sites. This PES is de- 
fined as (X,Y,Z,d,6, <f>)= {t/^{2),t/^{2),Z, 0.76A, 90°, 
135°), where t is the projected distance of the hydrogen 
center of mass from the sulfur atom in the (100) plane 
given in A. 

As one can see from Figure 6, the potential energy 
of the hydrogen molecule increases signficantly when 
the molecule comes close to the sulfur or the palladium 
atoms. The semi-circles in Figure 6 indicate the positions 
where the center-of-mass of the molecule is 1.5 A away 
from the S atom or 1.2 A away from the Pd atom, respec- 
tively. The potential energy of the hydrogen molecule at 
this distance has increased to values larger than 1.5 eV. 
This shows the strong inhibition of the hydrogen dissoci- 
ation in the vicinity of the sulfur and palladium atoms. 

Summarizing the results for the H2 dissociation on the 
(2x2) S covered Pd(100) surface, over the hi, h2, bi, 
tpd, tg sites we found that the dissociation behavior of 
H2 molecule strongly depends on the projected distance 
d|| between the hydrogen center of mass and the adsorbed 
sulfur atom, and on the orientation of the dissociating 
molecule. The dissociation behavior of hydrogen over the 
the (2x2) S covered Pd(100) surface can be summarized 
as follows: 

(1) For the geometry where the molecule reaches the 
surface at a position more than 3 A away from the sulfur 
atom , i.e. <i|| > 3 A (e.g. at the hi site), the dissociation 
is activated and the energy barrier is about 0.1 eV in the 
entrance channel resulting from the sulfur-induced mod- 
ification of the local electronic structure at the surface 
(see the elbow plots of bi-hi-bi and tpd-hi-tpd geome- 
try). 

(2) For the geometry where the molecule reaches the 
surface at a position about ~ 2.7 - 3.2 A away from the 
S atom (e.g. at the h2 or bi sites), if the molecule is 
oriented so that the hydrogen atoms do not approach 
the adsorbed sulfur, the dissociation is also activated. 
The energy barrier is in the entrance channel and its 
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magnitude depends on the distance of the molecule from 
the S atoms. For the hi-bi-li2 geometry, the energy 
barrier is 0.16 eV, whereas the barrier height is 0.6 eV 
for the b2-h2-b2 geometry because the latter adsorption 
position is closer to the sulfur atom. 

(3) For the geometry where the molecule reaches the 
surface at a position about ~ 2.7 - 3.2 A away from the 
sulfur atom (e.g. at the h.2 or bi sites), if the molecule 
is oriented so that the hydrogen atoms approach sulfur 
atoms, the dissociation becomes purely repulsive due to 
the direct interaction between H 2 and the S atoms (see 
the PES plots of b 2 -h 2 -b 2 and tpd-bi-tpd geometry). 

(4) For the geometry where the projected distance d|| is 
about ~ 1.8 - 2.3 A and if the molecule is oriented so that 
the hydrogen atoms do not approach the adsorbed sul- 
fur, the dissociation behavior of the hydrogen molecule is 
activated, however, with high energy barriers (> 1.2 eV) 
created by the strong repulsive interaction between hy- 
drogen and sulfur (see the PES of h2-tpd-h2 geometry). 

(5) For the geometry where the projected distance dii 
is smaller than 1.5 A, the dissociation behavior of the 
hydrogen molecule will be dominated by the strong re- 
pulsive interaction between H and S atoms and becomes 
activated with a high energy barrier (> 2.5 eV) (see PES 
of b2-ts-b 2 ) or purely repulsive if one of the H atoms 
points to the S atom. 

(6) For the geometry where the molecule is in an up- 
right geometry with its molecular axis perpendicular to 
the surface, the PES is repulsive. 

In Table [0] we summarize the results of the adsorption 
energies and dissociation barriers at six different geome- 
tries at the (2 x 2) sulfur covered surface and at the clean 
surface. For the clean Pd(100) surface, most of the hy- 
drogen dissociation pathways are non-activated if the H 2 
molecule dissociates parallel to the surface (0=90°), and 
only a few dissocation pathways are activated for this 
molecular orientation with small energy barriers (< 0.1^ 
eV). It was explained by Gross, Wilke, and SchefflcrEj 
that the large initial sticking coefficient (~0.7) oOpw 
energy hydrogen molecules on the Pd(100) surfaceE£l is 
mainly due to the steering effect. For the (2x2) sulfur 
covered surface, it is interesting to see that all the disso- 
ciation pathways become activated with energy barriers 
ranging from 0.1 eV to 2.55 eV. This fact explains qual- 
tatively the reason why the sticking coefficient of H2 is 
strongly reduced upon the adsorption of sulfur adatoms, 
however, a detailed uaderstanding needs a 6-D quantum 
dynamical calculationEa and will be published later .Ej 

In order to understand the different origin of the for- 
mation of small energy barriers at the hi and bi sites 
and large energy barriers at tpd, ts sites, we compare 
the density of states (DOS) for the H 2 molecule in dif- 
ferent geometries. Figure 7(a) shows the DOS when the 
H2 molecule is in the tpd-hi-tpd configuration. The H- 
H distance is taken as 0.75 A and the center of mass of 
the H 2 molecule is 4.03 A above the topmost Pd layer 
such that there is no interaction between the hydrogen 



molecule and the sulfur covered palladium surface. It is 
evident that the sulfur p orbitals strongly interact with 
the Pd d states, giving rise to a narrow peak just below 
the Pd d band edge (at e-ep — -4.8 eV) and a broad 
band at higher energies which has substantial DOS at 
the Fermi level. The d band at the surface Pd atoms is 
broadened due to the interaction with the S atoms. 

In Fig. 7(a), one intense peak of the DOS at the H 
atoms is found at the energy of the sulfur related bonding 
state at -4.8 eV. This degeneracy is accidental, however. 
There exists no broad distribution of states just below 
the Fermi level. This indicates that at this height above 
the surface, the a g orbital of the H 2 molecule interacts 
neither with the sulfur related bonding state nor with the 
surface Pd d band. 

When the H2 molecule gradually approaches the sur- 
face, the interaction of the hydrogen states with the sul- 
fur related bonding state and the surface Pd d states 
begins to build up. The intense peak of the DOS at - 
4.8 eV diminishs quickly and splits into a sharp bonding 
state and anti-bonding broad band just below the Fermi 
level. In Fig. 7(b) where the H 2 molecule is 1.61 A above 
the topmost Pd layer, the intense DOS of the H atoms 
found at -4.8 eV has shifted. The interaction of the a g 
orbital of the H2 molecule with the broad band of the 
surface Pd d states results in a sharp bonding state of 
the H2 CT 9 -surface interaction at -7.1 eV and a broader 
distribution of states with substantial weight below the 
Fermi level. Consequently we encounter an occupation 
of the H 2 -substrate antibonding states. Thus, a repul- 
sive contribution to the H2-surface interaction appears 
and gives rise to the formation of a small energy barrier. 

Figure 7(c) shows the DOS of the H 2 molecule in the 
b 2 -tg-b 2 geometry. The H-H distance is taken as 0.75 A 
and the center of mass of the H 2 molecule is 3.38 A above 
the topmost Pd layer. It is interesting to see that at this 
height a strong interaction is already found between the 
hydrogen molecule and the adsorbed sulfur. The intense 
peak of DOS at -4.8 eV has split into a sharp bonding 
state at (-6.6 eV) and a narrow anti-bonding state (at 
-4.0 eV) that strongly indicates the direct interaction of 
the H2 <T g orbital with the sulfur related 2 p state at -4.8 
eV which leads to a large energy barrier. 

Summarizing the above results and analysis, we 
conclude that the non-activated dissociation at clean 
Pd(100) surface is inhibited upon the adsorption of a 
(2x2) sulfur adlayer. The dissociative adsorption of hy- 
drogen with its molecular axis parallel to the surface is 
strongly corrugated; it has a wide range of energy barri- 
ers (0.10 - 0.15 eV near the hi and b2 sites; 0.6 eV near 
the h.2 site; 1.4 eV near the tpa site; 2.6 eV near ts) that 
strongly depends on the projected distance d\\ between 
H 2 and the S atoms, it also depends on the molecular ori- 
entation <j). Still dissociative adsorption of hydrogen is an 
exothermic process. Thus the poisoning effect of sulfur 
adatoms for H 2 dissociation at low sulfur coverage (83 < 
0.25) is governed by the formation of energy barriers and 
not by the blocking of adsorption sites. To our knowl- 
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edge, this is the most corrugated surface for dissociative 
adsorption studied so far by ab initio calculations. This 
has interesting consequences for the dissociation dynam- 
ics on this PES.El 



IV. H 2 AT THE c(2x2) SULFUR COVERED 
Pd(100) SURFACE 

A. Determination of the potential energy surface 
(PES) 

The potential energy surface for the dissociative ad- 
sorption of hydrogen over the sulfur covered c(2x2) 
Pd(100) surface is described within the same set of coor- 
dinates (A, Y, Z, dn-H, 6, (j)) as used for the (2x2) surface 
(see Figure 8). Figure 9 shows the surface unit cell for 
the c(2x2) sulfur covered Pd(100) surface. Inside the 
unit cell, we specifically analyse the dissociation at the 
following sites: bridge site (denoted by b), hollow site 
(denoted by h) , top site above the Pd atom (denoted by 
tpd), and top site above the S atom (denoted by ts). To 
obtain a detailed information about the interaction, wc 
have again calculated the elbow-plots at different sites 
and different molecular orientations (0,4>). 

B. Results 

For the c(2x2) S covered surface, the adsorbed S atoms 
form .a. 4 A x 4 A square lattice 1.29A above the Pd 
layer .til From that it follows that the projected distance 
d\\ of the molecule from one sulfur adatom cannot be 
larger than 2.85 A. Thus the hydrogen molecules on the 
Pd surface have to dissociate close to the sulfur atoms no 
matter how they approach the surface. Therefore, from 
what we have learned in the last section about the disso- 
ciation behavior of H2 over the (2x2) S covered surface, 
we might expect for the c(2x2) S covered surface that 
the dissociation will become purely repulsive for all dif- 
ferent sites and orientations. Indeed, this is what we find 
from our ab initio calculations. 



1. 2-D cuts through the PES at the hollow site 

At the h site, the projected distance d\\ between this 
site and the sulfur atom is 2.85 A. However, there are 
four S atoms near this site so that a dissociating hydro- 
gen molecule will always approach sulfur atoms. One can 
thus expect that the dissociation is purely repulsive for 
the hollow site. Figures 10(a), 10(b), and 10(c) show the 
PES for the t Pd -h-t Pd , the b-h-b, and the (X, Y, 8, 4>)= 
(0.5a, 0.5a, 0°, 0°) geometries at the hollow site, respec- 
tively. Figure 10 shows that the potential energy has 
increased to about 1.5 eV when the hydrogen molecule is 
at the same height above the Pd surface as the S atoms 



(.2 = 1.29 A). As mentioned above, this increase is due to 
the strong repulsive interaction between H and S atoms. 

2. 2-D cuts through the PES at top sites over Pd and S 
atoms 

At the tpd site, the projected distance d\\ to the nearest 
sulfur atom is only 2.0 A. Again we can expect that 
the strong repulsion between the H2 and the S atoms 
will dominate the dissociation behavior of the molecule. 
Figures 11(a) and 11(b) show the PES for the h-tp^-h 
and tg-tpd-ts geometries at the tp^ site. We find that 
the potential energy has increased to more than 2.0 eV 
even when the hydrogen molecule is at a height of Z=1.8 
A above the Pd atoms which is still 0.5 A above the sulfur 
adlayer. 

Figure 11(c) shows the PES at the ts site within the b- 
tg-b geometry. In this geometry, the hydrogen molecule 
approaches the surface directly towards the S atom. The 
potential energy increases dramatically to more than 1.5 
eV when the molecule is still 3 A above the topmost Pd 
layer. At this geometry the PES of the c(2x2) surface is 
found to be very similar to the PES of the (2x2) surface 
at the b2-tg-b2 geometry. We find an energy minimum of 
1.56 eV (at Z= 2.20 A, d H -H=l-90 A) as for the (2x2) 
surface. The fact that the PES of the (2x2) and the 
c(2x2) for this particular configuration are very similar 
shows that the hydrogen interaction with a sulfur covered 
surface close to the sulfur atoms is completely determined 
by the local bonding properties between the H and S 
atoms. 

Summarizing the results for the h, tp d , and tg sites, we 
find that the dissociation behavior of H2 at the c(2x2) S 
covered Pd(100) surface is dominated by the strong direct 
interaction between hydrogen and sulfur atoms. The PES 
for hydrogen dissociation is repulsive and endothermic 
due to the fact that for all different approach geometries 
the hydrogen molecules come close to the sulfur adatoms 
before the dissociation is completed. 

V. DEPENDENCE OF THE PES ON THE POLAR 
ANGLE 6 AND THE AZIMUTHAL ANGLE <j> OF 
THE HYDROGEN MOLECULE 

The two-dimensional cuts through the PES described 
in the section III and IV refer to the situation where the 
orientation of molecular axis is parallel or perpendicu- 
lar to the surface. In reality, the molecules impinges on 
the surface with all possible orientations. At the clean 
Pd(100) surface, Wilke and SchefneiQ have shown that at 
the h-b-h adsorption geometry the change of the orien- 
tation of the molecular axis with respect to the surface 
normal 9 away from 9 = 90° implies a significant increase 
in the potential energy. This increase was found n to be 
proportional to cos 2 9. Eichler, Kresse, and Hafncra have 
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found for the Rh(100) surface at less symmetric positions 
(halfway between bridge- and on-top position) that the 
potential energy also varies similiar to the cos 2 9 behav- 
ior but with the minimum of the energy shifted by 15° 
~ 20° from the parallel orientation of the molecule. 

To analyze the energy variation with 9 for the (2x2) 
sulfur covered Pd(100) surface, we first chose one high 
symmetry point (Z, c?h-h)= (1-05 A, 1.0 A) within the 
bi-hi-bi geometry and calculated the potential energy 
dependence on the angle 9 at this site. Figure 12 shows 
the energy variation with 9 which is well described by 
cos 2 9. To further analyze the energy variation with 9 
at less symmetric points, we chose the two points (Z, 
cZh-h) = (1.52A, 0.8A), and (Z, dn-n) = (1-05A, l.oA) 
at the hi-bi-h2 geometry and calculated their energy 
dependence on 9. The results are shown in Fig. 13. The 
potential energy is plotted as a function of cos 2 (6> — 9 a ) 
for values of 9 between 0° < 9 < 180°. Due to the lower 
symmetry, the potential energy is no longer single- valued 
in this representation as in the high-symmetry situation 
of Fig. 12 with 9 a — 0, but it is double- valued. The fact 
that the corresponding results fall upon one line indicates 
that the energy variation is well described by cos 2 (9 — 9q) 
with 9q ~ 2° - 5°. These findings are similar,-to the ones 
found for Rh(100) by Eichler and coworkersB 

The Figs. 12 and 13 demonstrate the significant in- 
crease in the potential energy by rotating the molecule 
away from 9= 90°, which is true at any point of the 
surface. This has been already indicated by the compar- 
ison of the PES of the hydrogen molecule approaching 
the hollow hi site in the tpd-hi-tpd geometry [see Fig. 
3(b)] with the PES of the hydrogen molecule at the same 
hi site in an upright position [see Fig. 3(d)]. It is evi- 
dent that rotating the molecule axis perpendicular to the 
surface increases the energy at the entrance channel and 
changes the dissociation pathway from being activated 
[see Fig. 3(b)] to being purely repulsive. 

The effect of the azimuthal angle <f> on the hydrogen 
dissociation can be seen by comparing the PES for the 
tpd-hi-tpd adsorption geometry [see Fig. 3(b)] with the 
bi-hi-bi adsorption geometry [see Fig. 3(a)]. For these 
two geometries at the hi site, the projected distance of 
the H2 center of mass from any sulfur atom is about 4 A 
and the two elbow plots are very similar at the entrance 
channel because the hydrogen atoms are still far away 
from any Pd or S atoms. However, differences are found 
at the exit channel because the hydrogen atoms disso- 
ciate towards the Pd atoms in the tpd-hi-tpd geometry 
whereas they dissociate towards the bridge site for the 
bi-hi-bi geometry. 

The dependence of the PES on the azimuthal angle 
is even more dramatic at adsorption sites close to the 
adsorbed sulfur. For example, if we compare the PES 
at the h 2 site in the b 2 -h 2 -b 2 geometry [see Fig. 3(c)] 
with the bi-h 2 -bi geometry [see Fig. 1(c) of Ref. 10], 
we can see that rotating the molecule towards the sulfur 
atoms leads to a large increase in the potential energy and 
changes the dissociation from being activated to being 



purely repulsive. In conclusion, we expect a significant 
increase of the potential energy by rotating the molecule 
towards sulfur atoms if the projected distance d\\ of the 
H 2 center of mass from the sulfur atoms is smaller than 
2.85A. 



VI. ANALYTICAL REPRESENTATION OF THE 

ab initio PES FOR H 2 AT S(2x2)/Pd(100) 

In order to perform dynamical calculation on the ab 
initio PES, one paeds a continous analytical representa- 
tion of the PES.E3 Details of such an representation for 
the PES of H 2 cm, the clean Pd(100) surface have already 
been publishedE3 For the H 2 dissociation on the (2x2) 
sulfur covered Pd(100) surface we have now also deter- 
mined an analytical representation which has been based 
upon the analytical PES for the clean surface. However, 
due to the larger unit cell of the S(2x2)/Pd(100) sur- 
face, the analytical form has been expanded up to the 
fourth-order Fourier cocffcient in the lateral directions 
parallel to the surface. In addition, in the azimuthal de- 
pendence of the PES the fourth-order term (proportional 
to cos(40)) has also been included. 

For the solution of the time-independent Schrodinger 
equation it is necessary to transform the coordinates, in 
the Zd-plane into reaction path coordinates s and p.EZI In 
these coordinates the function V(X, Y, s, p, 9, cf>), which 
describes the potential energy surface, has the following 
form in our parametrization: 

v (x,Y,s, P ,e,<f>) = v co " + v Iot + v vih (1) 

with 

4 

prcorr = ^ V$ n (s) cos mGX cos nGY, (2) 

m,n— 



V mt = V$ (s) - cos 2 9 (cos mGX + cos mGY) 

m— 

2 1 

+ Y1 V n 3 H s ) 2 sin2(9 cos2( ^ (cos nGX -cos nGY) 

n=l 

+ V (4) 0) \ sin 4 cos 40 (cos GX + cos GY) (3) 
and 

F vib = ^ w 2 (s) [p _ Ap{ x,Y,s)} 2 . (4) 

G = 2n/a is the length of the basis vectors of the square 
surface reciprocal lattice, a is the nearest neighbor dis- 
tant between S atoms in the (2x2)S/Pd(100) surface unit 
cell and uj(s) is the vibrational frequency. In addition, 
the curvature k = k(s) of the minimum energy path has 
to be determined. The displacement Ap in the potential 
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term V vlh (Eq. |^) takes into account that the location 
of the minimum energy path in the Zd-plane depends on 
the cut through the six-dimensional configuration space. 
Ap does not influence the barrier distribution, however, 
it changes the curvature of the minimum energy paths 
in the Zci-planes. Large values of Ap make the quan- 
tum dynamical calculations rather time-consuming since 
they require a large number of vibrational eigenfuactions 
in the expansion of the hydrogen wave functionO How- 
ever, the large values of Ap only occur for large separa- 
tions of the hydrogen atoms where they do not influence 
the calculated sticking probabilities and scattering prop- 
erties significantlyEH Therefore we have parametrized the 
displacement properly only for values of \Ap\ < 0.15 A. 

For the two hydrogen atoms in adjacent adsorption 
positions, the energetic cost to turn the molecule upright 
has not been determined by the ab initio calculations. 
This energy enters the term Vq in V rot (Eq. ||). Still 
it is possible to estimate that energy cost. In the har- 
monic approximation we can write the potential for the 
vibrations of the two hydrogen atoms perpendicular to 
the surface as 

Vperp0l,2 2 ) = yWad ( z l + 4)> ( 5 ) 

where z± and z-i denote the z-coordinates of the two 
hydrogen atoms, m is the mass of the hydrogen atom, 
and w a( j ist the vibrational frequency. Now we transform 
V'pcrp to the center-of-mass coordinate Z = {z\ + Z2)j2 
and the relative coordinate z — Z2 — Z\ . Then Vp Crp be- 
comes 

V peip (Z,z) = + f^ 2 (6) 

with the total mass M = 2m and the reduced mass p = 
m/2. The frequency w a d can be determined from the two- 
dimensional cuts of the ab initio calculations. Comparing 
V T0 t and Vp C rp and using z — cos (9 oIh-h, we obtain for 

(2) 

Vq at adjacent atomic adsorption positions 

Vq \s = Sad) = 2 W ad ^i?(ad)-if(ad) CD 

(2) 

Vq is therefore parametrized in such a way that it ap- 
proaches the value of Eq. ^ for the dissociated molecule 
on the surface. 

In general the functions Vm,n{s) and u>(s) are deter- 
mined such that the difference to the ab initio calcula- 
tions on the average is smaller than 50 meV. The energies 
are particularly well described along the minimum energy 
paths of the two-dimensional cuts of the six-dimensional 
configuration space. The maximum error of the analyt- 
ical representation in comparison to the ab initio PES 
is of the order of 0.5 eV. This error, however, does only 
occur for large distances from the minimum energy path 
where the potential energy is aleady rather high. There- 
fore the maximum error has only very little influence on 
the dynamical calculations. In Fig. Ilj we have plotted 



two cuts through the analytical six-dimensional PES for 
the bi-hi-bi and the h 2 -tpd-h 2 geometries, respectively. 
They should be compared with the corresponding ab ini- 
tio cuts in Fig. || and Fig. |[ 

VII. CONCLUSION 

We have performed detailed calculations of the PES 
for the dissociative adsorption of hydrogen molecules on 
the (2x2) and the c(2x2) sulfur covered Pd(100) sur- 
face using density-functional theory and-tlie full-potential 
linear- augmented plane wave method. BEJ The exchange 
correlation is treated in the generalized gradient approx- 
imation (GGA)B3 For the H2 dissociation on the (2x2) 
S covered Pd(100) surface, by calculating the PES over 
different adsorption sites we find that the non-activated 
dissociation pathways at the clean Pd(100) surface be- 
come activated or purely repulsive upon the adsorption of 
a (2x2) sulfur adlayer. The PES is strongly corrugated. 
The minimum barrier has a height of 0.1 eV, while close 
to the S atoms the barrier towards dissociative adsorp- 
tion for molecules with their axis parallel to the surface 
becomes larger than 2.5 eV. We find that the energy vari- 
ation with the polar angle 9 of the hydrogen molecule at 
the hi and bi sites is well described by cos 2 (8 — 6q) with 
6» « 0° - 5°. We also find that the PES strongly depends 
on the molecular azimuthal orientation <f> if the projected 
distance of H2 from the adsorbed sulfur is smaller than 
2.85 A. Still the dissociative adsorption of hydrogen is 
exothermic; thus the poisoning effect of sulfur adatoms 
for H2 dissociation at low sulfur coverage (0s < 0.25) 
is governed by the formation of energy barriers, not by 
blocking of the adsorption sites. 

For the c(2x2) S covered Pd(100) surface, the results 
of our calculations indicate that the interaction of H2 
with this surface is dominated by the strong repulsion 
between H2 and the sulfur atoms. On the c(2x2) S cov- 
ered Pd(100) surface the sulfur atoms form a 4 A x 4 A 
square lattice above the topmost Pd layer; therefore, all 
molecules that reach the surface will eventually approach 
some adsorbed sulfur atoms before the dissociation is 
complete. Due to this fact, non-activated reaction path- 
ways at the clean Pd(100) surface arc completely inhib- 
ited, and all dissociation pathways are purely repulsive. 
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FIG. 1. Sketch of a coordinate system for the descrip- 
tion of the hydrogen molecule above the (2x2) sulfur covered 
Pd(100) surface. The other coordinates are the height of the 
center-of-mass of H2 above the surface Z, the H-H distance 
dn-H, and the angle of the molecular axis with the surface 
normal 8. 



FIG. 2. Surface geometry of the (2x2) sulfur covered 
Pd(100) surface with two inequivalent hollow sites hi, \i2, 
two bridge sites bi, b2, top site tpd (above the Pd atom), 
and ts (above the S atom). 



FIG. 3. Cuts through the six-dimensional potential energy 
surface (PES) of H 2 dissociation over (2x2)S/Pd(100) at the 
hollow sites hi and I12: (a) PES for bi-hi-bi geometry; (b) 
PES for tpd-hi-tpd geometry; (c) PES for b2-h2-b2 geome- 
try; (d) PES for the molecule at the hi site with the molecular 
axis perpendicular to the surface. The energy contours, given 
in eV per molecule, are displayed as a function of the H-H 
distance, dij-H, and the height Z of the center-of-mass of H2 
above the topmost Pd layer. All length scales, also in the 
following, are given in A. The geometry of each dissociation 
pathway is indicated above. 



FIG. 4. Cuts through the six-dimensional potential en- 
ergy surface (PES) of H 2 dissociation over (2x2)S/Pd(100) 
at the bridge site bi: (a) PES for hi-bi-li2 geometry; (b) 
PES for tpd-bi-tpd geometry. The geometry of each dissoci- 
ation pathway is indicated above. 



FIG. 5. Cuts through the six-dimensional potential en- 
ergy surface (PES) of H 2 dissociation over (2x2)S/Pd(100) 
at the top sites tpd and ts: (a) PES for Ii2-tpd-h2 geome- 
try; (b) PES for b2-ts-b2 geometry. The geometry of each 
dissociation pathway is indicated above. 



FIG. 6. PES of a hydrogen molecule with d H -H= 0.76 A 
with its center of mass moving inside the (010) plane crossing 
the hi, tp d , ts sites. This PES is denned as (X,Y,Z,d,G, <f>) = 
{t/yf{2),t/yf{2),Z, 0.76A, 90°, 135°) with t and Z as the two 
variables given in A. The variable t is the projected distance 
of the hydrogen center of mass from the sulfur atom in the 
(100) plane. 



FIG. 7. Density of states (DOS) for a H2 molecule situated 
at (a) (Z,dif_ff)= (4.03A, 0.75A) within the tpd-hi-tpd ge- 
ometry; (b) (Z,d H -H)= (1.6lA, 0.75A) within the t Pd -hi-t P d 
geometry; (c) (Z,d H -ff)= (3.38A, 0.75A) within the b 2 -t s -b 2 
geometry. Given is the local DOS at the H atoms, the S 
adatom, the surface Pd atoms, and the bulk Pd atom. The 
energies are given in eV. 
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FIG. 8. Sketch of a coordinate system for the description 
of the hydrogen molecule above the c(2x2) sulfur covered 
Pd(100) surface. The other coordinates are the height of the 
center-of-mass of H2 above the surface Z, the H-H distance 
du-n, and the angle of the molecular axis with the surface 
normal 9. 

FIG. 9. Surface geometry of the c(2x2) sulfur covered 
Pd(100) surface with hollow site h, bridge site b, top site tp<j 
above the Pd atom, and top site ts above the S atom. 

FIG. 10. Cuts through the six-dimensional potential en- 
ergy surface (PES) of the H 2 dissociation at c(2x2)S/Pd(100) 
at the hollow site h: (a) PES for tp d -h-t Pd geometry; (b) PES 
for b2-h-b2 geometry; (c) PES for the molecule at the h site 
with its axis perpendicular to the surface. The geometry of 
each dissociation pathway is indicated above. 



FIG. 11. Cuts through the six-dimensional potential en- 
ergy surface (PES) of tch H 2 dissociation at c(2x2)S/Pd(100) 
at the top sites t Pd and ts: (a) PES for h-t Pd -h geometry; 
(b) PES for ts-ts-ts geometry; (c) PES for b-ts-b geome- 
try. The geometry of each dissociation pathway is indicated 
above. 

FIG. 12. Dependence of the potential energy of the H2 
molecule at the (2x2) sulfur covered Pd(100) surface on the 
angle 9 of its axis with the surface normal. The potential 
energy is plotted as a function of cos 2 9. The configuration for 
9 — 90° corresponds to the bi-hi-bi geometry with (cIh-h, 

z) = (1.0A, 1.05A). 

FIG. 13. Dependence of the potential energy of the H2 
molecule at the (2x2) sulfur covered Pd(100) surface on the 
angle 6 of its axis with the surface normal shifted by 8 . The 
potential energy is plotted as a function of cos 2 (9 — 8 ) for 
values of 6 between 0° < 9 < 180°, i.e. the potential energy 
is double-valued in this representation. The configuration for 
9= 90° corresponds to the hi-b2-h2 geometry: (a) (cIh-h, Z, 
9 ) = (1.52A, 0.8A, 4.5°); (b) {d H -H, Z, 9 a ) = (l.oA, 1.05A, 
2.2°). 

FIG. 14. Two cuts through the analytical representation of 
the six-dimensional ab initio PES for the bi-hi-bi and the 
h 2 -tp d -h2 geometries, respectively. 
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geometry 


(Z, d) 


3 Pd slab, £ cut = 11 Ry 


5 Pd slab, E cut = 14 Ry 


bi-hi-bi 


(0.61 A,0.97 A) 


0.28 eV 


0.24 eV 


hi-bi-h2 


(0.08 A,2.85 A) 


-0.79 eV 


-0.72 eV 


(0.97 A,1.14 A) 


0.13 eV 


0.08 eV 


tpd-hi-tpd 


(1.05 A,0.94 A) 


-0.02 eV 


-0.01 eV 


h2-tp d -h 2 


(1.45 A, 1.29 A) 


1.30 eV 


1.31 eV 



TABLE I. Dependence of the potential energy on the thickness of the substrate Pd slab and the cut-off energy at different 
H 2 geometries. The energies are given per H2 molecule. 



adsorption site and geometry 


clean surface (0 S =O) 


(2x2) surface (6 S =0.25) 


hi (0,0)=(9O O ,9O°) 


non-activated, E ad = 0.46 eV 


activated E b = 0.10 eV 


h 2 ((9,0)=(9O O ,9O°) 


activated E b = 0.60 eV 


h 2 (0^)=(9O°, 0°) 


repulsive (> 0.75 eV) 


hi (<9,<rt=(90 o ,135°) 


non-activated, E a d = 0.23 eV 


activated Eb= 0.13 eV 


hi (0,0)=( 0°, 0°) 


no dissociation 


repulsive (> 0.25 eV) 


bi (e,0)=(9O°,9O°) 


non- activated, E a d= 1.22 eV 


activated Eb= 0.15 eV 


b 2 (9,<t>)=(90°, 0°) 


no dissociation 


repulsive (> 0.25 eV) 


t Pd (6»,0)=(9O°,135°) 


activated, E b = 0.16 eV 


activated E b = 1.28 eV 


t s (6,^(90°, 0°) 




activated E b = 2.55 eV 



TABLE II. Summary of the results for adsorption energies and barrier heights for the H 2 dissociation over the clean Pd(100) 
surface and the (2x2) sulfur covered Pd(100) surface at different adsorption sites and geometries. The energies are given per 
H 2 molecule. 



11 




Fig. 2 






Fig. 5 




Fig. 7 




Fig. 9 




Fig. 11 





Fig. 13 




Fig. 14 



